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We show that to n loop order the divergent content of a Feynman amphtude is spanned by a set 
of basic (logarithmically divergent) integrals Jjog(A^), i = 1, 2, • ■ ■, n, A being the renormalization 
group scale, which need not be evaluated. Only the coefficients of the basic divergent integrals are 
show to determine renormalization group functions. Relations between these coefficients of different 
loop orders are derived. 
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I. INTRODUCTION 



Implicit regularization (IR) is a non-dimensional momentum space framework which has been claimed 
to be a strong candidate for an invariant regularization suitable to develop perturbation theory in super- 
symmetric gauge field theories []|-(T7|. Assuming an implicit regulator in a general (multiloop) Feynman 
amplitude, a mathematical identity at the level of propagators allows to write the divergent content as 
basic divergent integrals (BDI) or loop integrals written in terms of one internal momentum only in an 
unitarity preserving fashion. This is possible because BPHZ subtractions as well as the counterterm 
method are compatible with IR to arbitrary loop order. An arbitrary scale appears via a regularization 
independent identity which relates two logarithmically BDF s by trading a mass parameter m (or an 
infrared regulator in the propagators) for an arbitrary positive parameter A, ([A] = M) plus a function of 
to/ A. Consequently A parametrizes the freedom of separating the divergent content of an amplitude and 
acts as a renormalization group scale. The key point underlying IR is that neither the (regularization de- 
pendent) BDF s nor their derivatives with respect to A represented by BDFs need be evaluated. In other 
words, the BDFs are readily absorbed into renormalization constants whose derivatives with respect A 
used to calculate renormalization group functions can also be expressed by BDF s. The advantage of such 
scheme is that a physical amplitude is written as a finite part plus a set of BDFs say //^^(A^) and finite 
surface terms (ST's) expressed by volume integrals of a total derivative in momentum space which stem 
from (finite) differences between //^^(A^) and /«^i^=^ -(a2) where the latter is a logarithmically divergent 
integral which contains in the integrand a product of internal momenta carrying Lorentz indices /ii, /i2---- 
In other words throughout the reduction of the amplitude to loop integrals, ^/og'^^'^^ "(A^) may be written 

as a product of metric tensors symmetrized in the Forentz indices times I^lgi^'^) plus a surface term. 

Such ST's are in principle arbitrarily valued. However it has been shown that setting them to zero ab 
initio corresponds to both invoking translational invariance of Green's functions and allowing shifts in the 
integration variable in momentum space [E] which in turn is an essential ingredient to demonstrate 
gauge invariance based on a diagrammatic proof. Therefore ST's seem to encode the possible symmetry 
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breakings. Moreover it has been verified that constraining such surface terms to nought is also sufficient 
to guarantee that supersymmetry is preserved in tlie Wess-Zumino model to 3'"^-loop order and 
supergravity to 1-loop order [3| • Notwithstanding it is reasonable to assert that IR is a good candidate 
to an invariant calculational friendly regularization framework valid in arbitrary loop order. From the 
point of view of algebraic renormalization, ST's would be the necessary symmetry restoring counterterms 
whose expression is known within IR. Then a constrained version of IR (CIR) amounts to setting them 
to zero from the start and thus constituting an invariant scheme. When physical quantum breakings 
(anomalies) are expected some care must be exercised: one is able to spot a genuine breaking by letting 
the ST's to be arbitrary so to verify that none consistent set of values for the ST's dictated by symmetry 
requirements fulfill all the essential Ward identities of the underlying model at the same time 0] , ■ In 
, 01 the rules that define IR to arbitrary loop order are specified. 

A renormalization group equation can immediately be written within IR adopting A as a renormal- 
ization group scale and a minimal, mass independent renormalization scheme in which only the basic 
divergent integrals are absorbed in the renormalization constants. Hence renormalized Green's function 
satisfy a kind of Callan-Symanzik equation governed by the scale A. 

The purpose of this contribution is to twofold. Firstly although IR works in arbitrary massive quan- 
tum field theories, for massless theories it undergoes a remarkable simplification. Assuming an infrared 
regulator fi for the propagators, llggip"^) equals -f;og(A^), (A ^ 0), plus a sum of terms proportional to 
powers of the logarithm of the ratio fJ,/X. We will show in this contribution that for massless theories all 
the divergencies to arbitrary loop order can be cast as a function of ^/^^(A^), according to the definition 

where = J (d'^k)/ {2Tr)^ and the superscript A is a symbol for an implicit regularization. Secondly 
it is well known that renormalization group functions constitute a testing ground for regularizations 
because they both encode the symmetry properties of the underlying model which should be preserved 
by the regularizations and their expansion in perturbation theory contains terms which are universal, i. 
e. renormalization scheme independent. While some interesting simplifications take place in dimensional 
methods, e.g. in an inverse power series in e — > of the coupling constant, beta functions are determined 
uniquely by the residue of the simple pole on e, it is pertinent to ask what is the counterpart in IR. 
That is to say, one may wonder how the calculation of renormalization group functions systematizes 
within a scheme where only basic divergent integrals are claimed to be suHicient to exhibit the ultraviolet 
properties of a model in a symmetry preserving fashion. The answer to this question is that a general 
framework for renormalization group functions can be built in which the simplifications of dimensional 
methods manifest themselves as relations between the coefficients of basic divergent integrals coming 
from different Feynman graphs that contribute to a given renormalization group function. 

We illustrate with the Yukawa model in 3 -I- 1-dimensions to 2"''-loop order which contains a 75 matrix 
and hence the application of dimensional regularization is more involved. 



II. GENERAL ULTRAVIOLET STRUCTURE OF MASSLESS THEORIES 

The purpose of this section is to show that the ultraviolet content of an amplitude to n*^ loop order 
for massless models, considering the definition, is written in terms of llogi^"^)- A general n-loop, Z-point 
amplitude, after space-time and internal group algebra contractions are performed, can always be written 
as a combination of integrals of the type 



(fc • • • (fc -pi)' 
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where we have integrated n — 1 times leaving only fc, the most external loop momentum and the p^'s are 
external momenta. For a massless model suppose that An-i is cast like 

" / fc2\ - 

i—l ^ ^ 

in which An-i is finite under integration on k and A^-i, the divergent part, represents the subdivergences 

which in principle are already written in terms of ^/^^(A^). The mass scale has emerged from a 
scale relations which characterizes a renormalization scheme in Implicit Regularization. The coefficients 
ai{k,pi, ■ ■ ■ ,pi) may contain powers in the external and internal momenta. To justify the assumption of 
equation ^ we proceed with a proof by induction. For n = 2 (one loop order) it can be easily verified 
that ([31) holds for Ai . Now we show that this assumption for (n — l)*''-loop order implies the same 
structure for the n*''-loop order to conclude by induction that the multiloop integrals at any order have 
the same structure. The relevant contributions come from the second term on the r.h.s. of ([3|), 



/' 

Jk 



[{k-p^Y~^?]■■■[{k~plY-^i 
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ln"M-^^T^ . (4) 



which has superficial degree of divergence r{i) — 2/ + 4. Extra factors in the numerator were considered 
so as to account for the Lorentz structure of the ai{k,pi, ■ ■ ■ ,p;)'s. A fictitious mass fj,^ was introduced 
in the propagators and the limit /^^ — > will be taken in the end. A fictitious mass may always be 
introduced if the integral is infrared safe. This is necessary because although the integral is infrared 
safe, the expansion of the integrand, as we explain below, breaks into infrared divergent pieces. When 
a genuine infrared divergence appears, this procedure can be problematic in non-abelian theories. For 
such cases a new procedure within IR defining basic infrared divergent integrals is necessary in order to 
preserve symmetries [l3'|. 

We judiciously apply in the integrand the identity. 



1 1 pi.-2pr-k 



[Pr - ky - fl^ (fc2 - /i2) (fc2 _ ^2) [(p^ _ /,)2 _ ^2] 

for the factor in the denominator which depends on pi to obtain 



(5) 



r(2)-2/+5 



m=l 



V (-1)'"-^ / — ' ' ' In , 

^ K (fc2-^2)m[(^_p^)2_^2]...[(fc_pj_^)2_^2] \^ ^2 



^ ^ A (fc2-^2).W-2/+5[(fc_p^)2_^2]...[(^_p^)2_^2] ^2 ^ " ^ ^ 

In the equation above the last integral is finite. We pick out the ultraviolet divergent ones. For each 
one of the divergent integrals in the summation, the procedure has to be repeated for all the external 
momenta. Let us consider one typical divergent integral after the expansion has been performed for all 
the external momentum with the exception of p\ : 

■J.^ -,r-J^ (fc2_^2)a[(fc_p^)2„^2]il^ ^2 )■ 

For this integral the superficial degree of divergence is r — 2a + 2 and the expansion is performed so as 
to have the divergent part freed from the external momentum: 



m=l 



(fc2 - ^2)a+m \^ X 



\2 



^ ' Jk (fc2-^2).-o+3p_p^)2_^2] 1^ ^2 )■ 
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For the basic divergent integrals (without dependence on the external momenta), it is only possible to 
have an even degree of divergence. Besides, as shown in ref. (from equations (20) to (23) of this 
reference), it is possible to write a parametrization in which the quadratic divergences vanish in the limit 
/i^ — > to one loop order. The same argument can be generalized to arbitrary loop order. So, we only 
have to deal with the basic logarithmic divergent integrals. They have the form, 



^ fc'^i • • • k^'^ 

(fc2 - ^2)p 



In* 



A2 



(9) 



where p = a + m, m being the summation index in equation ([H]), and 2p = j + 4, which in turn may 

always be written in terms of -?^;og(A*^)'s (see equation ([T])) plus surface terms. For example, for two 
Lorentz indices we have 



+2(i 




A2 



(10) 



The procedure is repeated for I^^^ ^'^^^ so to obtain 



1 O'-l)! 
2J-' {i - l)!"'°f 



IillifJ-'') + surface terms. 



(11) 



We still have to deal with the fictitious mass, which in the limit ^^2 ^ will give infrared divergent 
pieces both in the ultraviolet divergent and finite parts. This problem is simply dealt with by the use of 
regularization independent scale relations (they can be easily obtained with the help of a cutoff), which 
read 



■E 

fe=i 



(J-I)! 
A:! 



In" 



(12) 



for arbitrary non- vanishing A. This justifies the appearance of the mass scale A2 in An-i- For infrared 
safe models a systematic cancelation of all powers of In ^ between the ultraviolet divergent and finite 
parts finally crowns A a renormalization group scale. The important fact here is that an integral of the 
type ([7]) will have a general result given by 



In™ ( 



(13) 



where T^™.'.^^ and are tensor structures depending on the external momentum. This permits us 

to conclude that in the next loop order, the same structure of divergence will be maintained. The other 
finite parts, obtained when the expansion of the integrand was carried out, when inserted into an external 
loop can give divergent contributions. Nevertheless, they can be put in the form of ([7]) by performing 
further expansions. Therefore we can assert that the divergent structure of massless loop calculations 
within the context of Implicit Regularization can be completely displayed in terms of the /;^*^(A2)'s. This 
completes our proof. 

For the sake of clarity, we exemplify below: 



1 



/c2(fc-p)- 



■In'^ 



2p- fc 



, (fc2-Ai2)2p_p)2_^ 




A2 



(14) 
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FIG. 1: Feynman rules to the Yukawa model: (a) ^ ij (b) i/p^, (c) —ig, (d) —675, (e) ^ iAp^ , 

(f) *B (g) ^ -iOff, (h) ^ -675^ 



which for n = 3 yields 



hm * y-2! 

2^0 I 167r2 fc! 



fc=U 

where in the last step we have used (IT^ for j = 3 





i 




167r2 


fc! 


/ p2 


\ A2 



A;=0 



P2 

A2 



(15) 



III. RENORMALIZATION GROUP FUNCTIONS 



In this section we present a general framework to work out renormalization group functions using 
renormalization constants defined by BDI's. We will see that derivatives of BDI's which are also BDI's 
need not be evaluated. Moreover, in the calculation of renormalization group functions, a simplification 
becomes manifest through relations between some coefficients of BDI's. We study the massless Yukawa 
theory in 3 + 1-dimensions to 2-loop order as a working example because both it is rich enough due to 
the presence of overlapping divergences and two coupling constants. Besides, dimensional methods are 
more involved as a 75 matrix appears in the interaction term. 

11 1 
The Lagrangian density in terms of renormalized variables (pQ = Z^cj), ipo — Z^^tp, cq = eZe/{Z-^Z^), 

90 = gZg/Z"^ reads 

£ = (1 + A)^^(|)^(|)^' + i{l + B)4,-i^d^'TP + i{l + C)e7^75i/;0 - (1 + (16) 

where Z^ = l + A, Z.^ = l + B,Ze = l + CcZg = l + D. 

The superficial degree of divergence for any graph with external boson lines and external fermion 
lines is given by A = 4 — — thus to 2 loop order the divergent amplitudes are {n^ = 2 and 
= 0), {n^ = and = 2 ), T^^^^ (n^ = 1 e n^, = 2 ) and r^4 {n^ = 4 and = ) . Next we 
evaluate the diagrams necessary to compute the renormalization group functions to 2-loop order, which 
are portraited in figures 2 to 5. It is not difficult to show that the amplitudes following from the Feynman 
rules can be treated in IR by separating the external momentum dependence in the BDI's using ([5|) and 
neglecting surface terms which stem from ()lip whilst the renormalization group scale is defined through 
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FIG. 2: Diagrams contributing to V ^i . 
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FIG. 3: Diagrams contributing to F 



(fT^ . We summarize the results below. In figure 2 the divergences can be isolated as 



2c 



r2b 
r2d 

r2f 

r2h 



2eV/io3(A2) 



4 2 

e p 

4 2 

e p 
167r2 

4 2 r 

e p 



167r2 

4 2 

e p 



87r2 

2 2 

5 P 



-5IiogiX^) + zl6^2 [/;^^(a2)]% 2/z„g(A2)ln f-^ 



12(47r) 



and r 



02 



iAp^ + r 



2b + 2c 
2 



2d 



167r2 



2g + which yields 



2c + -L 2f 
^2 \ 



(17) 
(18) 

(19) 

(20) 

(21) 

(22) 

(23) 
(24) 

(25) 



In determining the counterterm graphs that correspond to the amplitudes in equations (|19p and (|20p 
above we have used the one loop contributions of the counterterms B and C from equation ([29]) . Notice 
that the non-local divergences have been correctly canceled as they should because we have shown in 
that IR is compatible with the counterterm method derived from BPHZ forest formula. Diagrams in 
figures 3, 4 and 5 are evaluated in a similar fashion, using one loop counterterms previously determined, 
to give, 



167r2 



(26) 



■ TpTp(j) 



167r2 



(27) 



r^4 = -igD + {2AiT:'^g'^ - 3847r2e^ - 6e^ - 125^62 + 336e^ + 9656*^ 



167r2 



+ (3ff^ + 65^6^ - 144e6 - 7256*) + , / _ _ gg^e _ ^^^^a | [j,^^(a2)]2^ (33) 
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FIG. 4: Diagrams contributing to V 
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FIG. 5: Diagrams contributing to Tj^i. 



respectively. 

Then we have the renormaUzation constants defined in a minimal scheme if 



A 

B 

C 
D 



167r2 

i 

167r2 
1 



12 



31 



167r2 



e2l67r2/,„,(A2) + eV;o,(A^) + (9//o,(A^) - 6/g(A')) 
^b247r^/,„,(A2) +5^ (6z7;o,(A2) - Si/gCA^) - 12n-'[ho,{\^)f) 
+ g-ie4384i7r2/,<,3(A2) + .g-ie*^' (-SSGzI/oglA^) + 144z4'](A2) - Q&{l&TT^)[Ii^g{X^)] 
+ {-miio^iX') + 12il[ll{X^) - 72(167r2)[/,„,(A2)]2)], 
The beta-functions and field anomalous dimensions are defined as usual 



A^ dZ^ _ dZ^ 

^* ~ z^'d^'^^~'z^~d^' 

/3e = -ffA^ (Z-^^-^Z 



dx^ 2^ dx^ ^' aA^j'^s- ^^-^ [^s 



-idZ^X 



(29) 



(30) 



where A is the IR arbitrary scale which plays the role of renormalization group scale. To n-loop order, a 
general renormalization constant can be written as 



(31) 



in which p, r, s assume positive integer values in each n. In a minimal, mass independent renormalization 



scheme, the renormalization constants Z^^ , Z^^ and z'^g} take the general form 

Z« = ^ A«[/,o,(A^)]'= + E4^'^4'j(A^)- 



(32) 
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For the Yukawa model to two loop order we have, 

2 



z„ 



l + g^Zi^) + Y.^'-Z, 

Tt = ] 

1 + ^c"^^'"' 

n=l 

n=l 

1 + e^zf + ge^Zf + J] (g-ie^^+^zi") + 5"^, 



o2n 



(n) 



(33) 



Now plugging equations (j33p into pop permits us to obtain the finite contributions to renormalization 
group functions to 1 and 2-loop order from 



(1) _ „2\2^^^6^ „,(!)_ „2\2C^ 



(1) 



— e A 



2^3,2/1^ dzP dzl\ 

\2 9A2 dX^ dX^ )' 

p,r,{l) f,y{l) 07(1) 



(34) 



and 



(2) 



2"^a ^p4"-b ) (2)^g4;^2^^ 



.2e35A2^ + eg^A^^ + 2e'X^ ' ''^^ 



1 dzP ^ 



9A2 



aA2 



9A2 2 9A2 9A2 



aA2 aA2 



2g^A 



3 \2 



dzi^^ dz. 



(2)~ 



aA2 aA2 



2e^X'^^'^ 



(2) 



aA2 



2 2 2x2£^;s_ 



(35) 



To complete our task we have to evaluate the derivatives of ([32]) w.r.t A2 which are expressible in terms 
of BDFs as well, namely 



A" 



^dZ^ 



dX^ 



167r2 



j=2 



fc=2 



(36) 

a = a, ■ • ■ I, which is a general expression for massless models though a similar one holds for massive 
models as well. Direct inspection of equation enables us to determine the coefficients yl^"-* and i?^"-* 
which appear in (j34p and ([35]). After some straightforward algebra we obtain 



1 



3 



48- 



2 (47r)2 ' (47r)2 ' (47r)2 (47r)2 (47r)2 



(37) 
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and 

= (2)_ 13 

"^"^ 12 (47r)4 (47r)4 ' ^V' g (47r)4 ' 

" "T (4^ + 12 (4^ " (4^' ^» " (4^ + (4^ ~ (4^ + ^(4^' ^'^^^ 

which agree with [l^, p^ . 

We can generahze ([M)) and ([55]) to arbitrary loop order using the expansion (pij) in ([50)1 to conclude 
that all we need to evaluate the renormalization functions is the derivative of Z*^"^ as given in It 
is interesting to remark that whilst the finite terms in the r.h.s of (|36p contribute to the computation 



of the renormalization group functions, the terms proportional to BDI's will give relations between A^"'' 

and S^""* as they must vanish because the renormalizations functions are finite. The same reasoning 
leads us to conclude, in dimensional regularization methods, that only residues of order one contribute 
to beta-functions. For instance, from the calculation of the field anomalous dimensions 7^ and 7^ up to 
two loop order we get. 



A^a^-^^^^B^;^^ = (39) 



( 



4? - m-B^S) = \ (<■)' + 4i' {-\a'S + 4!' - A'S ] . (40) 



and 



. (a^ - S^.^B^^) ^ \ {a^^)\ + A« - 41)) (41) 

respectively. 

To conclude we have shown that in Implicit Regularization (IR) , we can organize the divergent content 
of an amplitude to n*'* loop order in terms of a basis of basic divergent integrals (BDI's), namely {Iiog{^^)} , 
i = 1 ■ ■ - n where A is the RG scale. The calculation of RG functions systematizes within IR for they can 
be written in terms of coefficients of BDI's. Such coefficients are shown to be inter-related which in turn 
allows us to restrict ourselves to a subset of BDI's at each loop order to evaluate RG functions. 
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Appendix 

We calculate explicit the diagram g of figure 2: 

r2, = -try^ i (-■'^^^)/-^^^)i(-^^^)^(--'^^'^)F^(fc37)I (42) 
where I and k are internal momenta. Taking the trace of Dirac matrices and simplifing we obtain: 



r2, = -2z5^ / / 

Jk Jl 



I Pk'^(k~pY{i~pY{k-i)^ 



(43) 
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or 

^■2,--2^9'^^~p'l k^i^k-prl w^^^L k^l {i-pnk-iy} 

At this point we apply in each of these integrals the methods discussed in section II. After some algebra 
we get: 

Tag = ^ (^-Miog{X')+il67r^ [liogiX")]^ + 2IiogiX'')lii (^-^^ + finite^ (45) 

Observe that the third term on the r.h.s. of P5)) is non-local and it must be canceled with the ones of 
the countcrtcrm diagrams. 
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